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Abstract 

In this paper we study the finite time blow-up problem for the ax- 
isymmetric 3D incompressible Euler equations with swirl. The evolu- 
tion equations for the deformation tensor and the vorticity are reduced 
considerably in this case. Under the assumption of local minima for 
the pressure on the axis of symmetry with respect to the radial varia- 
tions we show that the solution blows-up in finite time. If we further 
assume that the second radial derivative vanishes on the axis, then 
system reduces to the form of Constantin-Lax-Majda equations, and 
can be integrated explicitly. 
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1 The axisymmetric 3D Euler equations 

We are concerned with the following Euler equations for the homogeneous 
incompressible fluid flows in a domain Q C M 3 , 

div v = 0, (1.2) 



1 



v(x,0) = v Q (x), (1.3) 
where D/Dt is the material derivative defined by 

Here v = (vi,V2,v 3 ), Vj = Vj(x,t), j = 1,2,3, is the velocity of the flow, 
p = p(x, t) is the scalar pressure, and vq is the given initial velocity, satisfying 
div Vq = 0. Since the classical result ([12]) on the local well-posedness for the 
3D Euler equations in the standard Sobolev space H m (M. 3 ), m > 5/2, the 
problem of finite time singularity for such local smooth solution is still an 
outstanding open problem, although there is a celebrated result on the blow- 
up criterion ([1 J) and its refinements (pD, M, E]) taking into account geometric 
considerations on the vorticity directions. By an axisymmetric solution of the 
Euler equations we mean velocity field v(r, x 3 , t), solving the Euler equations, 
and having the representation 

v(r, x 3 , t) = v r (r, x 3 , t)e r + v e (r, x 3 , t)e e + v 3 (r, x 3 , t)e 3 

in the cylindrical coordinate system, where 



e r ={— ,— ,0), e e = {- — ,— ,0), e 3 = (0,0,1), r=Jx\ + x 2 2 . 

In this case also the question of finite time blow-up of solution is wide 
open(see e.g. [21 El H] for preliminary studies of the problem; see also [10] for 
the related recent result in case of helical symmetry) . The vorticity u = curl 
v is computed as 



where 



lu = tu r e r + uj e eo + u> 3 e 3 , 



u r = -d 3 v 9 , J = d 3 v r - d r v 3 , u 3 = — + d r v e . 



r 

We denote 

v = v r e r + v 3 e 3 . 
The Euler equations for the axisymmetric solution are 

(v e ) 2 



r 



d t v r + (v ■ \7)v r = —d r p + 
d t v 3 + (v ■ V)v 3 = -d 3 p, 

r 8 

d t v e + (v-V)v e = , 

r 

div v — 0, 
k v(r,x 3 ,0) = v Q (r,x 3 ), 

where V = e r d r + e 3 d 3 . Note that the above representation of the Euler equa- 
tions in the cylindrical coordinate system is valid off the axis of symmetry, 
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which is chosen to be x 3 — axis. Hence, in order to analyze the equation 
on the £3— axis, we mainly use the equations in the Cartesian coordinate 
system. Below the functional values of a cylindrically symmetric function 
f(xi, X2, x 3 ) = f(r, x 3 ), r = \fx\ + x\ on the £3— axis, should be understood 
as Hm r _ 40 + f{r, x 3 ) := f{x 3 ). 

Theorem 1.1 In the axisymmetric 3D Euler equations with the symmetry 
of axis chosen to be x 3 —axis, we write 

u ■= d x v 2 - d 2 v 1 (= 2d lV 2 = 2d r v 9 ), A := d 3 v 3 , 

which are defined on the x 3 —axis. Suppose the initial data satisfies 

S := {x 3 e R I u> (x 3 ) = 0, \ {x 3 ) > 0, d 2 rPo {x 3 ) > 0} ^ 0, 

where we denoted 

u} (x 3 ) = u)(x 3 ,0), A (x 3 ) = A(x 3 ,0), p (x 3 ) = p(x 3 ,0). 

We define T\ = Ti(a) as 

T x = inf{t > I d?p(X 3 (a,t),t) < 0}, 



where X 3 (a,t) is the particle trajectory defined by the local classical solution 
v(x, t). 

0X:i{( '- f) v 3 (X 3 (a,t),t), X 3 (a,0) = a. 



dt 

Then, there exists no global classical solution to the axisymmetric 3D Euler 
equations if there exists a G S such that 

T\(a) > -ly (1.4) 
A (a) 

Remark 1.1 After the above theorem is proved P. Constantin informed me 
of the preprint ( [6 J ) , where it is shown that the positivity of all the matrix 
components of the hessian of the pressure leads to a singularity in the general 
case. In the above, however, the positivity is assumed essentially only for 
one component on the x 3 — axis in the case of axisymmetry. 



Remark 1.2 The assumption on the positivity of the second radial deriva- 
tive of the pressure on the x 3 — axis is physically natural in view of the fol- 
lowing heuristic argument. We consider a axisymmetric compressible ideal 
fluid with swirl. Due to centrifugal force the density of fluid becomes local 
minimum on the x 3 — axis, which implies local minimum of pressure on the 
x 3 — axis, hence d 2 p > on the axis. Now we take zero Mach number limit for 
the pressure to obtain the pressure of the original axisymmetric incompress- 
ible fluid (see [13 J for rigorous result for this singular limit problem). In this 
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limiting procedure it is plausible to expect preservation of the local minimum 
property ofthe pressure on the x 3 — axis. 



Theorem 1.2 In the axisymmetric 3D Euler equations with the symmetry 
of axis chosen to be x 3 — axis, let us assume that there exists T > such that 

d 2 r p{x 3l t) = V(x 3 , t)eRx [0, T]. (1.5) 

Then, the pair (u, X), which is defined in Theorem 1.1, can be explicitly given 
by 

(2 - X {a)t) z + u Q [a)H 
\(Y In A A 4Ao(q)-2[Ao(a) 2 + ^o(a) 2 ]t 
(2 - X {a)tj z + u {a)H 

along the particle trajectory {X 3 (a,t)} for all (a,t) G R x [0, T]. Let us 
assume S = {x 3 G | Ao^) > 0, ^0(0:3) = 0} 7^ 0. Then, the form of 
solution fli.7] ) implies that there exists no global classical solution to the 3D 
axisymmetric Euler equations if 

T > inf } 



aes A (a) 

2 Proof of the main theorems 

We begin with the following elementary lemma. 

Lemma 2.1 Let v = (v ,v ,v 3 ) = v(x\,x 2 ,x 3 ) be an axially symmetric 
C 1 — vector field a on R 3 with the axis of symmetry chosen as the x 3 — axis, sat- 
isfying div v = 0, and let p = p{x\, X2, x 3 ) be an axially symmetric C 2 — scalar 
function on R 3 . Then, on the axis of symmetry we have 

v 1 = v 2 = 83V 1 = d 3 v 2 = d lV 3 = d 2 v 3 = 0, (2.1) 
v r = v e = d 3 v r = d 3 v 9 = d r v 3 = 0, (2.2) 



d lV l = d 2 v 2 = = d r v r = lim — , (2.3) 

2 r^o r 

3 lV 2 = -82V 1 = d T v e = lim — , (2.4) 

r- >o r 

dip = d 2 p = did 2 p = did 3 p = d 2 d 3 p = d r p = d r d 3 p = 0, (2.5) 

d\p = d 2 2 p = clip = lim (2.6) 
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Proof Here we use notations, 



T T / 

(x 1: x 2 ) := (-x 2 ,xi) and (x l: x 2 ) := (—^=, -y=), r=^x\ + x 2 2 . 

Let us observe first 

v 1 (xi,x 2 ,x 3 ) = — v r - — v 9 , v 2 (x 1 ,x 2 ,x 3 ) = — v r + — v 9 , r > 0, 

ry ny ry ry 

and thus 

v 1 (x[,x' 2 ,x 3 ) = -—v r - yf 9 = -v 2 (x 1 ,x 2 ,x 3 ), (2.7) 

v 2 {x[, x 2 , x 3 ) = —v r -v 9 = v 1 (x\, X2, x 3 ). (2.8) 

Passing r — > in (I2.7p - fl2.8l) . we find that V\ — v 2 — on the x 3 — axis. On 
the other hand, 

v 1 (x 1 ,x 2 ,x 3 ) +v 2 (x 1 ,x 2 ,x 3 ) = V2v r (xi,x 2l x 3 ), (2.9) 
v 1 (xi,x 2 ,x 3 ) - v 2 (xi,x 2 ,x 3 ) = -V2v e (x 1} x 2 ,x 3 ) } (2.10) 

and passing r — > in (j2.9j) - (12.10j) . we also find that v 9 = v r = on the 
x 3 — axis. Replacing v\, v 2 , v r , v 9 by d^ 1 , d 3 v 2 , d 3 v r , d 3 v 9 respectively in the 
above argument we also deduce that d 3 v 1 = d 3 v 2 = d 3 v r = d 3 v 9 = 0. Next 
we note that 

div 3 (xi,x 2 ,x 3 ) = ^d r v 3 , d 2 v 3 (xi,x 2 ,x 3 ) = ^-d r v 3 , r > 0, 



and therefore 



div 3 (x' 1 ,x' 2 ,x 3 ) = -—d r v 3 = -d 2 v 3 (xi,x 2 ,x 3 ), (2.11) 

r 

d 2 v 3 (x[,x 2 ,x 3 ) = —d r v 3 = d 1 v 3 (x 1 ,x 2 ,x 3 ) (2.12) 

r 

for all r > 0. Similarly to the above, passing r — > in fl2.11l) - fl2.12l) . we 
deduce d\V 3 = d 2 v 3 = on the x 3 — axis. Since 

div 3 (x 1 , x 2 , x 3 ) - d 2 v 3 (xi,x 2 , x 3 ) = V2d r v 3 (x 1 , x 2 , x 3 ), 

we are lead to d r v 3 = on the 23— axis by passing r — > 0. In order to verify 
(12.31) we compute 

r 2 2 

q 1/ \ ^ ^1 r 1 x l a r , ^l 2 ^ ^1^2 # 

div (xi,^,^) = nV + ^rd r v H ^-ti 5-a r w , 

ry ryO ry£ ryrJ ry-di 

r 2 2 

O 2/ \ ^ ^2 r , ^2 q r ^1^2 , a 'l a '2 

(9 2 f (xi,x 2 ,x 3 ) = -v + — d r v —v H —d T v 

iy ryC) ty £i ryO ry £i 
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for r > 0. Hence, 



div (x[,x' 2 ,x 3 ) = d 2 v (x 1 ,x 2 ,x 3 ) Vr > 0. 



(2.13) 



Passing r — > in (I2.13p . we have d\v l = d 2 v 2 on the x 3 — axis. The condition 
div v = implies dif 1 = c^f 2 = — ^c^f 3 . We note 

<9i?/(zi, x 2 , x 3 ) - <9 2 v 2 (xi, x 2 , x 3 ) = ! — ! d r v 9 (x 1 ,x 2 , x 3 ), (2.14) 

r 

and 

^it) 1 ^!,^,^) + d 2 v 2 (x 1 ,x 2 ,x 3 ) = ! — ! h (9 r i7 r (xi,x 2 ,a;3). (2.15) 



lim — = lim d r v 9 . 

r— »0 r r— >0 



(2.16) 



$2 2 

^1^2 r , X\X 2 v x 2 e x 2 „ e 



2 2 
XiX 2 Q r f X 1 g X 1 e 



nr> ryttj 



From (I2.14p we have 

Let us compute 

d 2 v 1 (x 1 ,x 2 ,x 3 ) = 
d 1 v 2 (x 1 ,x 2 ,x 3 ) = 
and, hence 

div 2 (x[,x' 2 ,x 3 ) = -d 2 v 1 (x 1 ,x 2 ,x 3 ) (2.17) 

for all r > 0. Passing r — > in (12.171) . we obtain <9it> 2 = — d 2 v x on the 
X3— axis. Let us compute 

d 2 V 1 (xi,X 2 ,X 3 ) +diV 2 (xi,X 2 ,X 3 ) = — - ( Xl ' X2 ' X3 ) +g r t;' r (x 1 ,X2,X3), (2.18) 

r 

and 

d 2 v 1 (x u X2,x 3 ) -div 2 (xi,X2,x 3 ) = — d r v e (x 1 ,x 2 ,x 3 ). (2.19) 



The equation (12.181) provides us with 



lim — = lim d r v r , 

i — >0 T r^O 



(2.20) 



while (12.191) . combined with (12.161) . shows div 2 = d r v e respectively on the 
23— axis. Using the fact (12.201) . passing r — > in (12. 15ft . we deuce d\Vi = d r v r 
on the x 3 — axis. As for (12.51) the proof of dip = d 2 p = d\d 3 p = d 2 d 3 p = d r p = 
d r d 3 p = is exactly same as the above and we omit it. We note 



d 1 d 2 p{x 1 ,x 2 ,x 3 ) 



X\X 2 XiX 2 2 

— —d r p+ —^-d r p, 



and find that did 2 p(x[, x' 2 , x 3 ) = — did 2 p(xi, X2, x 3 ). Hence, passing r — > 0, 
we have <9i<9 2 p = on the x 3 — axis. We also compute 

d' 2 p(x 1 ,x 2 , x 3 ) = -d r p - -^d r p + ~^d 2 r p = d\p{x x , x' 2 , x 3 ), 

and deduce that 

d\p = d\p on the x 3 — axis by passing r — > (2.21) 

Note that 

d\p + d 2 2 p = -d r p + d 2 r p, (2.22) 
r 

and 

X 1 

dxd 2 p(xi, x 2 , x 3 ) = ~i^prP + ^rP as r ^ 0. (2.23) 
From fl2T2TD . (12T22D and (12T231 we have d\p = dip = d 2 .p = lim^o \d r p. □ 

Next we recall the matrix representation of the Euler equations(see e.g. |14j). 
Given velocity v(x,t), and pressure p(x,t), we introduce the 3x3 matrices, 

v _ dvj_ V i:j + Vjj _ Vjj - Vjj d 2 p 



ctei 2 2 dxidxj 

with z,j = 1, 2, 3. Then, we have the decomposition V = (Vy) = S+A, where 
S = (Sij) represents the deformation tensor of the fluid, and A = (Aij) is 
related to the vorticity u by the formula, 



^3 3 
A^ = - '^^EijkUk, &i = ZijkAjk, (2.24) 
fe=i j,k=i 



where e^k is the skewsymmetric tensor with the normalization e 123 = 1. Note 
that P = (Pij) is the hessian of the pressure. Let {Ai, A 2 , A 3 } be the set of 
eigenvalues of S. Computing partial derivatives d/dxk of (11.11) yields 

^— = -V 2 - P. (2.25) 
Dt y ' 

Taking symmetric part of (12.251) . we have 

D ' S ' -S 2 -A 2 -P, (2.26) 



Dt 

from which, using the formula ( 12.24(1 . we derive 
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DSii _ 

Dt • 4 

fe=i 



^ SikSkj + -{\uj\ 2 5ij - LOiUJj) - P^, (2.27) 



where 5^ is the Kronecker delta defined by 5^ = 1 if i = j, and 5^ = 
otherwise. The antisymmetric part of (12 .251) is 

DA 

— = -SA - AS, (2.28) 
Dt v ' 

which, using the formula (12.241) again, we obtain easily 

lkC Su, (2.29) 



Dt 



which is the well-known vorticity evolution equation that could be derived 
also by taking curl of (11.11) . Taking trace of (12.271) . we have the identity 

- (A? + Al + A^) + i|^| 2 = Ap. (2.30) 

Proof of Theorem 1.1 Thanks to lemma 2.1 we have the following reduced 
representation for the deformation tensor, the vorticity, and the hessian of 
the pressure on the x 3 — axis. 

S = diag(-~,-^,A), u = (0,0, a;), P = di ag (d 2 p, d 2 rP , d 2 3 p), (2.31) 
where 

A = d 3 v 3 (= ~d lV l = -U 2 v 2 = -\d r v^j , 00 = d lV 2 - d 2 v l (= 2d r v e ) 

on the 23— axis. The (11) and (22) components of the matrix equation (12.271) 
reduce to 

DX X 2 lu 2 o2 . 

m = T - T + 2d ' p ' (2 ' 32) 

where we set 

§- t = d t + v 3 d 3 , 
while the (33) component becomes 

^ = -A 2 - dip. (2.33) 
We note that f!2.30j) reduces to 

Ap = -^A 2 + ^, (2.34) 

which is also obtained by taking subtraction (I2.32I) - (I2.33I) . The vorticity 
equation is written as 

^ = AO, (2.35) 
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which can be solved as 



uj(X 3 (a, t),t) = u> (a) exp 



X(X 3 (a, s), s)ds 



along the trajectory. This implies that uj(X 3 (a,t),t) = for a G S as long 
as classical solution persists. Hence, for a G S H2.33j) can be written as 



d\(X 3 (a,t),t) X 2 (X 3 (a,t),t) 



dt 



> 



X 2 (X 3 (a,t),t) 



+ 2%p(X 3 (a,t),t) 

Vt G (0,11(0)). (2.36) 



The differential inequality (12.361) can be solved immediately to yield 
X(X 3 (a,t),t) > 2X °^ Vt G (0,7;) with T» =T*(a) := min{Ti(a ~ 



2-tAnfa) 



A (a) 



which shows that Ti(a) > 2/Ao(a) is not consistent with the fact that classi- 
cal solution persists until T(a). □ 



Proof of Theorem 1.2 By the hypothesis the equation f!2.32j) together 
with (I2.35P reduces to 

DX _ A 2 _u 2 _ 
Dt 2 2 ' 

Duj _ 

= Xuj. 

Dt 

This is exactly the same system studied by Constantin-Lax-Majda in [S] 
with the material derivative replacing the partial derivative in time, which 
was proposed as a one dimensional model equation for the 3D Euler equations 
in the vorticity formulation. Similarly to [8] we set = A + iu. Then (12.371) 
becomes the following complex Riccati equation along the trajectory, 

De _ e 2 

~Dt ~ ~2~' 

which can be solved explicitly as 

e(X 3 (a, t),t) = = g^#±^L . (2 .38) 

v 6\ , h ) 2-e {a)t 2 - [A (a) + iu Q (a)]t K ' 

Taking imaginary and real parts of (12.38}) we obtain (ll.6p -( |L7l) . □ 

Remark after the proof In [UJ Hou-Li also obtained a system of equa- 
tions similar in form to (12.371) . but for a different pair of unknown functions 
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under completely different assumptions. In our case the system is derived rig- 
orously from the axisymmetric 3D Euler equation by taking the limit r — > 0, 
and assuming only dip = on the X3— axis. 
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